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1 Introduction 



The horizon and flatness problems of the standard big-bang cosmology are solved el- 
egantly if during the evolution of the early universe the energy density happened to 
be dominated by some form of vacuum energy which resulted in a quasi-exponential 
growth of the scale factor [|TJ . An inflationary stage is also required to dilute any unde- 
sirable remnants such as topological defects surviving from some phase transition at a 
pre-inflation epoch. 

The vacuum energy driving inflation is generally assumed to be associated with some 
scalar field 0, known as the inflaton, which is initially displaced from the minimum of 
its potential. As a by-product of solving the horizon and flatness problem, quantum 
fluctuations of the inflaton field may produce the seeds necessary for the generation of 
structure formation and for fluctuations in the cosmic background radiation. 

Inflation ended when the potential energy associated with the inflaton field became 
smaller than the kinetic energy of the field. By that time, any pre-inflation entropy 
in the universe had been inflated away, and the energy of the universe was entirely 
in the form of coherent oscillations of the inflaton condensate around the minimum of 
its potential. The universe may be said to be frozen after the end of inflation. We 
know that somehow the low-entropy cold universe dominated by the energy of coherent 
motion of the </> field must be transformed into a high-entropy hot universe dominated 
by radiation. The process by which the energy of the inflaton field is transferred from 
the inflaton field to radiation has been dubbed reheating.^] 

In the old theory of reheating |2|], it was assumed that the inflaton field oscillated 

around the minimum of its potential in a coherent way until the age of the universe 

grew to the order of the inflaton decay lifetime r^, t ~ = FT . At this stage, 

the inflaton decayed and the universe filled with the inflaton decay products, which 

soon thermalized. In the process the universe was "reheated" to the temperature of 

Tjih — lO^yT^Mpi, where Mpj ~ 10 19 GeV is the Planck mass. In a simple chaotic 

inflation model the inflaton potential is given by V(<p) = m 2 2 /2, with m ~ 10 13 GeV in 

1 Reheating may well be a misnomer since there is no guarantee that the universe was hot before in- 
flation. Since we are confident that the universe was frozen at the end of inflation, perhaps "defrosting" 
is a better description of the process of converting inflaton coherent energy into entropy. 
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order to reproduce the observed temperature anisotropies in the microwave background. 
Writing = a^m, one finds T RH ^ lO 15 ^^ GeV [§]. 

The density and temperature fluctuations observed in the present universe, Sp/p ~ 
10 -5 , require the inflaton potential to be extremely flat. This means that the couplings 
of the inflaton field to the other degrees of freedom (including a^) cannot be too large, 
since large couplings would induce large loop corrections to the inflaton potential, spoil- 
ing its flatness. As a result, Trh is expected to be smaller than 10 14 GeV by several 
orders of magnitude. The problem of large loop corrections to the inflaton potential 
may be solved in the framework of supersymmetry H, where the nonrenormalization 
theorem || guarantees that the superpotential is not renormalized at any order of per- 
turbation theory. On the other hand, in supergravity-inspired scenarios gravitinos have 
a mass of order a TeV and a decay lifetime on the order of 10 5 s. If gravitinos would 
be overproduced in reheating and decay after the epoch of nucleosynthesis, they would 
modify the successful predictions of big-bang nucleosynthesis. This can be avoided if 
the reheat temperature is smaller than about 10 11 GeV (or even less, depending on the 
gravitino mass). 

In addition to entropy, the baryon asymmetry must be created after inflation. One 
method to generate the baryon asymmetry is by the decay of baryon-number (B) violat- 
ing superheavy bosons (referred to generically as "X" bosons, whether gauge or Higgs 
bosons) of Grand Unified Theories (GUT's) |J. In the old theory of reheating there 
is a serious obstacle to post-inflation GUT baryogenesis related to the relatively large 
X-boson mass and the relatively small reheat temperature. 

The unification scale is generally assumed to be around 10 16 GeV, and 5-violating 
gauge bosons should have masses comparable to this scale. Baryon-number violating 
Higgs bosons may have a mass one or two orders of magnitude less. For example, in 
SU (5) there are B violating "Higgs" bosons in the five-dimensional representation that 
may have a mass as small as 10 14 GeV. In fact, these Higgs bosons are more likely than 
gauge bosons to produce a baryon asymmetry since it is easier to arrange the requisite 
CP violation in the Higgs decay f7|, |8|, [J. But even the light 5-violating Higgs bosons 
are expected to have masses larger than the inflaton mass, and it would be kinematically 
impossible to create them directly in decay. 
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One might think that the X bosons could be created by thermal scattering during 
the stage of thermalization of the decay products of the inflaton field. However, if Trh 
is as small as necessary to avoid overproduction of gravitinos, production of superheavy 
bosons by thermal scattering would be heavily suppressed.f] 

But the outlook for GUT baryogenesis has brightened recently with the realization 



that reheating may differ significantly from the simple picture described above [\ \ 



13| , |14| , 15 1 . In the first stage of reheating, called "preheating" [fLlJ| , nonlinear quantum 
effects may lead to an extremely effective dissipational dynamics and explosive particle 
production even when single particle decay is kinematically forbidden. Particles can 
be produced in the regime of a broad parametric resonance, and it is possible that a 
significant fraction of the energy stored in the form of coherent inflaton oscillations 
at the end of inflation is released after only a dozen or so oscillation periods of the 
inflaton. What is most relevant for the present discussion is that preheating may play 
an extremely important role for GUT generation of the baryon asymmetry. Indeed, it 
was shown in |16|] that the baryon asymmetry can be produced efficiently just after the 
preheating era, thus solving many of the problems that GUT baryogenesis had to face 
in the old picture of reheating. 

A crucial observation for baryogenesis is that even particles with mass larger than 
that of the inflaton may be produced during preheating. To see how this might work, 
let us assume that the interaction term between the superheavy bosons and the inflaton 
field is of the type g 2 (p 2 \X\ 2 . During preheating, quantum fluctuations of the X field with 
momentum k approximately obey the Mathieu equation: X£ + [A(k) — 2qcos2z]Xk = 0, 
where q = g 2 (p 2 /4m 2 , A(k) = (k 2 + m 2 x )/m 2 + 2q, and primes denotes differentiation 
with respect to z = int. Particle production occurs above the line A = 2q. The width of 
the instability strip scales as q 1 ' 2 for large q, independent of the X mass. The condition 
for broad resonance, A — 2q < q 1 ! 2 |TTJ] , becomes (k 2 + m 2 x )/m 2 < g(j)/m, which yields 
for the typical energy of X bosons produced in preheating E\ = k 2 + m 2 x < gcf)m [T3|] , 
where <f) is the amplitude of the oscillating inflaton field. By the time the resonance 



2 There exists another problem for GUT baryogenesis scenarios: B violation through sphaleron 
transitions are expected to be fast at high temperatures, and would erase any preexisting baryon 
asymmetry produced at the GUT scale |lfj unless there is a non vanishing value of B — L. But a 
natural way to overcome this problem is to adopt a GUT like 5*0(10), where an asymmetry in B — L 
may be generated. 
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develops to the full strength, 2 ~ 10~ 5 Mp 1 . The resulting estimate for the typical 
energy of particles at the end of the broad resonance regime for m ~ 10 _6 Mpi is Ex ~ 
10 _1 g 1 / 2 i/mMpi ~ (7 I / 2 10 15 GeV. Supermassive X bosons can be produced by the broad 
parametric resonance for Ex > %, which leads to the estimate that X production will 
be possible if m x < # 1/2 10 15 GeV. 

For g 2 ~ 1 one would have copious production of X particles (in this regime the 
problem is non-linear from the beginning and therefore g 2 = 1 has to be understood as a 
rough estimate of the limiting case) as heavy as 10 15 GeV, i.e., 100 times greater than the 
inflaton mass.[] The only problem here is that for large coupling g, radiative corrections 
to the effective potential of the inflaton field may modify its shape at <ft ~ -^pi- However, 
this problem does not appear if the flatness of the inflaton potential is protected by 
supersymmetry. 

This is a significant departure from the old constraints of reheating. Production of 
X bosons in the old reheating picture was kinematically forbidden if m < mx, while 
in the new scenario it is possible because of coherent effects. It is also important to 
note that the particles are produced out-of-equilibrium, thus satisfying one of the basic 
requirements to produce the baryon asymmetry [|l^ . 

Scattering of X fluctuations off the zero mode of the inflaton field limits the maximum 
magnitude of X fluctuations to be (X 2 ) max ps m 2 / g 2 |14]]. For example, (X 2 ) max ~ 
10 _10 M|j in the case mx = 10 m. This restricts the corresponding number density of 
created X-particles. 

A potentially important dynamical effect is that the parametric resonance is effi- 
cient only if the self-interaction couplings of the superheavy particles are not too large. 
Indeed, a self-interaction term of the type A|X| 4 provides a non-thermal mass to the X 
boson of the order of (A(X 2 )) 1//2 , but this contribution is smaller than the bare mass 
mx, if A < g 2 m 2 x /m 2 . Self-interactions may also terminate the resonance effect because 
scattering induced by the coupling A may remove particles from the resonance shells 



and redistribute their momenta [16], [19| . But this only happens if, again, A g [13 



The parametric resonance is also rendered less efficient when the X particles have 
a large decay width Y X - Roughly speaking, one expects that the explosive production 



3 In the case in which the cross-coupling between the inflaton and the X field is negative, superheavy 
particles may be produced even more efficiently |]l7| . 
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of particles takes place only if the typical time, r e , during which the number of X 
bosons grows by a factor of e, is smaller than the decay lifetime tx = T^ 1 . During the 
broad resonance regime, typically r e > 10 m -1 . If we parameterize the decay width by 
Tx = amxi this requires a < 0.1m /mx 0- Notice that smaller values of Tx are favored 
not only because particle production is made easier, but also because the superheavy 
particles may remain out-of-equilibrium for longer times, thus enhancing the final baryon 
asymmetry. 

The exact knowledge of the maximum allowed value of the decay width of the su- 
perheavy degrees of freedom is therefore of extreme importance for the computation of 
the final baryon asymmetry produced by the GUT particles after preheating. 

The goal of this paper is twofold. First, we wish to provide the first fully non- 
linear calculation of the inflaton decay into superheavy X bosons taking into account 
their decay width Tx- Our basic finding is that the parametric resonance does not 
develop if the decay rate Tx is larger than about 10 _1 m, thus confirming the rough 
estimate made above. X production through the resonance is very efficient for smaller 
values of Tx- Our second goal is to compute numerically the number density nx of 
supermassive X bosons produced at the resonance stage. This parameter is fundamental 
for the computation of the final baryon asymmetry. We will also show that as long as 
the bound Tx ^ 10 _2 m is satisfied, the observed baryon asymmetry B ~ 4 x 10~ n 
may be explained by the phenomenon of GUT baryogenesis after preheating, with no 
further restriction of the parameters. We will also comment on the phenomenological 
implications of our findings. 



2 X production and decay 



Using the methods developed in Refs. |T2], [TB|, [TJ|], we have studied numerically the 
production of massive, unstable X particles in the process of the inflation decay. 

We consider a model in which the oscillating inflaton field <fi interacts with a scalar 
field X whose decays violate baryon number B. A simple possibility for the X-particle 
is the Higgs field in the five- dimensional representation of SU(5), although as noted 
above, because of the desirability of B — L violation 5*0(10) is a more promising theory. 



^eealsoref. |0| 



5 



We assume standard kinetic terms, minimal coupling with gravity, and a very simple 
potential for the fields of the form 

V(<P, X) = ImV + \m\X* + 1/0 2 X 2 . (1) 

Let us first introduce dimensionless variables in the conformal reference frame. 
The rescaled conformal time r is related to cosmological time t by mdt = a(r)dr. 
The rescaled conformal fields x an d f are related to the original fields by X{r) = 
x(r)0 o (O)/a(r) and <f>(r) = <^(r)0 o (O)/a(r). In this model <j> (0) w 0.28M P1 is the 
value of the inflaton field at the end of inflation. We assume that immediately af- 
ter inflation the universe is matter dominated and the scale factor evolves as a(r) = 
(V^>o(0)t/V3M p1 + l) 2 @. 

An important dimensionless parameter in the problem will be the resonance param- 
eter q = 3 2 2 (O)/4m 2 . For m = 1.3 X l(r 6 M P1 and </> (0) = 0.28M Pb q ~ 10 10 ^ 2 . 

In the conformal variables, the equations of motion become 

(p - VV + (a 2 -^)<P+ 4 «V = , 
X - V 2 x + Tax + (rn 2 x a 2 -Td-^j X + = . (2) 

Here we are taking into account the decay of the X field by simply introducing the term 
Tx X in the equation of motion for the field X. The dimensionless parameter T which 
enters the Eq. (0) is T = Tx / m. Similarly, m x = mx / m. 

We have solved these equations of motion directly in coordinate space on a 128 3 
spatial lattice with periodic boundary conditions. The initial conditions for fluctuations 
correspond to the conformal vacuum at the time when the oscillations of tpo commence 
[T2"| , |T3|| . The initial conditions for the coherently oscillating inflaton zero- momentum 
mode are y?o(0) = 1, 0o(O) = 0. 

A fundamental parameter in GUT baryogenesis is nx, the number density of the 
supermassive leptoquarks whose decays produce the baryon asymmetry. It will depend 
upon the value of T and q. 

Since the supermassive bosons are more massive than the inflaton, one expects small 
kinetic energy in the excitations of the X field. From the potential of Eq. (0), the square 
of the effective mass of the X field is (m x FF ) 2 = (m x + g 2 (<f) 2 )) and the energy density 
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Figure 1: The variance of X with model parameters q = 10 6 , m x = 2, and T = 6 x 10 -2 is 
shown by the lower solid curve as a function of time. The upper solid curve corresponds 
to the inflaton zero mode. The dotted curves represent the same quantities for T = 0. 

in the X field will be p x = {m 2 x + g 2 (4> 2 )) (X 2 ) . Writing (0 2 ) as 0q + (<50 2 ), we can 
define an analog of the X-particle number density as 



1/2 



n x = Px/mT = [Hr + W 2 »/0o(O) + m 2 j m(X 2 } . (3) 

Eq. (0) enables us to calculate the number density of the created X-particles if the 
variances of the fields, (X 2 ), (5<j) 2 ) = ([(/>(x) — (4>)} 2 ), and the inflaton zero mode 0o( r ) = 
(0(r)), are known. Here (■ • •) has to be understood as the average over statistical 
realizations. Since the system is homogeneous on average, this is equivalent to the 
volume average. We shall present the dependence upon the time of the variances for 
several choices of the model parameters as well as the maximum value of the variances 
which may be achieved during the evolution as a function of the same parameters. 

The time evolution of the variance, (X 2 ), and of the inflaton zero mode, (0), is 
shown in Fig. [TJ by the solid curves for the case q = 10 6 , m x = 2, and T = 6 x 10~ 2 . 
We see that the particle creation reaches a maximum at r w 10.8 when (X 2 ) « 10 -9 
in the "valleys" between the peaks.P| At later times, r > 10.8, particle creation by 



3 Note that we can use Eq. (g) only for the "valley" values of the variance, where the adiabatic 
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Figure 2: The maximum value of the variance of the X-field, (X^ ax ), is shown as 
function of T. Stars mark (X^ ax ) obtained in the full non-linear problem. (X^ ax ) in 
the Hartree approximation is shown by the dotted curve for q = 10 4 , m x = 0.1, and by 
the solid curve for q = 10 6 , m x = 2. 



the oscillating inflaton field can no longer compete with X-decays due to the non-zero 
value of T. For comparison, we show in the same figure the case T = represented by 



the dotted curves [fLl}] . In the T = case, particle creation is able to compete with the 
expansion of the universe so that (X 2 ) remains roughly constant. Another novel feature 
of the case with non-zero T is the low level of inflaton fluctuations, (5<J) 2 )(t) <C <Pq{t). 
Using Eq. @, we find for the maximum number density of created X-particles 



nx 



4g^(lO.8)/0g(O) + m? m(X 2 ) 



1/2 



10 3 + ml 



1/2 



m(X 2 ) w 30m(X 2 ). 



It is easy to understand that if we increase the value of T, the parametric resonance 
will not be able to compete with the decay of X at earlier times. Moreover, for suffi- 
ciently large values of T, the resonance will be shut off in the linear regime. One goal 
of this paper is to find the boundary of the model parameter space that will result in 
sufficient X-particle creation for successful baryogenesis. 

In exploration of parameter space we turn to the Hartree approximation (for details 



approximation is valid, while non-adiabatic amplification occurs in the region of the peaks of (X 2 (t)) 
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see Ref. |13|), which requires much less computing resources. The maximum value 
of the variance of X reached during the time evolution of the fields in the Hartree 
approximation is shown in Fig. ^| as a function of the parameters of the model. Here the 
stars also show the maximum of (X 2 (t)) in the full non-linear problem for a few values 
of T. At small T the Hartree approximation overestimates (X 2 ) significantly [fl3| , |14}| . 
Nonetheless, at large values of T it is a quite reliable approach. We see that (X 2 ) drops 
sharply when T > 0.2, and we have checked that this critical value of T does not depend 
significantly upon mx or q. 

The most relevant case with q = 10 8 , where X-bosons as massive as ten times the 
inflaton mass can be created, is shown in Fig. 3 in the Hartree approximation. Note, 
that two lower curves which correspond to T equal to 0.08 and 0.12 never reach the 
limiting value (X 2 ) max ~ 10 ^-/kfpj, which is imposed by rescattering ||14|| , and the 
Hartree approximation ought to be reliable in this cases. 

The final baryon asymmetry depends linearly upon the ratio 5 between the energy 
stored in the X particles at the end of the preheating stage and the energy stored in 



the inflaton field at the beginning of the preheating era fl6fl . From our results, we can 
estimate that this ratio as 

^3xl0^ m ,<g. (4) 

Therefore, for q = 10 8 and m x = 10, 5 is of the order of 3 x 10 8 (X 2 )/M^ 1 . Since 
the final baryon asymmetry scales approximately as T _1 and is given by B ~ 5 x 
10 _4 5e(r/5 x 10~ 5 ) -1 [|T(J, where e is an overall parameter accounting for CP violation 
(it will be typically a one-loop factor times some CP- violating phases), we see that the 
observed baryon asymmetry B ~ 4 x 10~ n may be explained by the phenomenon of 
GUT baryogenesis after preheating if 



M 2 X \ e ) V5 x 10- 5 

From Fig. 3 we can read that this only may happen if Tx is smaller than about 10~ 3 mx- 
This result may be considered very comfortable since we can conclude that whenever 
the resonance develops, i.e., when Tx ^ 10 _1 m = 10 _2 mx, GUT baryogenesis after 
preheating is so efficient that the right amount of baryon asymmetry is produced for 
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T 

Figure 3: The time dependence of the variance of X in the Hartree approximation with 
model parameters q = 10 8 , m x = 10 and for three values of T, from top to bottom: 
0.04, 0.08, 0.12. 

almost the entire range of values of the decay rate Fx- In other words, provided that 
superheavy X-bosons are produced during the preheating stage, they will be ineffec- 
tive in producing the baryon asymmetry only if their decay rate falls in the range 
10- 3 m x <T X < 10- 2 m x . 



3 Conclusions 

From our findings we may infer some phenomenological implications for any model 
of GUT baryogenesis in preheating. As we already mentioned, B violation through 



sphaleron transitions are expected to be fast at high temperatures pi , and would erase 
the baryon asymmetry produced after preheating unless the supermassive X bosons 
generate some nonvanishing value of B — L. A natural way to overcome this problem 
is to adopt a GUT like 5*0(10), where an asymmetry in B — L may be generated by 
the lepton- number violating decays of the Higgs field 0i 2 6 in the 126-representation of 
SO (10) which transforms as a singlet under the SU(5) decomposition. 
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This Higgs field is responsible for the Majorana masses of the right-handed neutrinos 
through the symmetric couplings h UR (pi26 • VrCvr. Requiring that T^, 126 is smaller than 
about lCr 3 m^ 126 imposes the constraint h UR < 10 _1 . Through the see-saw mechanism, 
light neutrinos may therefore get masses in the range 10~ 3 to 10~ 2 eV for h VR ~ 10~ 3 . 

Notice that, as opposed to large- angle scattering processes, forward-scatterings do 
not alter the distribution functions of the particles traversing a gas of quanta, but only 
alter the dispersion relation. This remains true even in the case of a nonequilibrium 
system such as the one represented by the gas of </>i 2 6 particles created in preheat- 
ing. At the end of reheating the right-handed neutrinos receive a mass of the order of 

(0126(1)) 1 I 2 from non-thermal corrections with the 0i26(l)-particles due to forward- 
scatterings. However, this mass is not large enough to suppress kinematically the de- 
cay rate of the 0i26- This would require the energy of the 0i26's to be smaller than 
h VR (0i26(l)) 1 / 2 , which only happens if ((0i26(l))/-^Pi) Xfi larger than about 10~ 8 g h~£. 

In conclusion, we have performed a numerical analysis of the production of super- 
heavy X-bosons during the preheating stage following the end of inflation and have 
shown that the observed baryon asymmetry may be produced in the decay of these 
non-thermal GUT bosons if the value of their decay rate is smaller than about I0~ 3 mx- 
GUT baryogenesis after preheating solves many of the serious drawbacks of GUT baryo- 
genesis in the old theory of reheating where the production of superheavy states after 
inflation was kinematically impossible. Moreover, the out-of-equilibrium condition is 
naturally attained in our scenario since the distribution function of the X-quanta gen- 
erated at the resonance is far from a thermal distribution. This situation is considerably 
different from the one present in the GUT thermal scenario where superheavy particles 
usually decouple from the thermal bath when still relativistic and then decay producing 
the baryon asymmetry. 

It is quite intriguing that out of all possible ways the parametric resonance may de- 
velop, Nature might have chosen only those ways without instantaneous thermalization 
and also with a successful baryogenesis scenario. 
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